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Abstract 


Three control law design techniques for flutter suppression are pre- 
sented. Each technique uses multiple control surfaces and/or sensors. The 
first method uses traditional tools (such as pole/zero loci and Nyquist dia- 
grams) for producing a controller that has minimal complexity and which 
is sufficiently robust to handle plant uncertainty . The second procedure 
uses linear combinations of several accelerometer signals and dynamic 
compensation to synthesize the modal rate of the critical mode for feed- 
back to the distributed control surfaces . The third technique starts with a 
minimum- energy linear quadratic Gaussian controller , iteratively modifies 
intensity matrices corresponding to input and output noise, and applies 
controller order reduction to achieve a low-order , robust controller. The 
resulting designs have been implemented digitally and tested subsonically 
on the active flexible wing wind-tunnel model in the Langley Transonic Dy- 
namics Tunnel. Only the traditional pole/zero loci design was sufficiently 
robust to errors in the nominal plant to successfully suppress flutter dur- 
ing the test. The traditional pole/zero loci design provided simultaneous 
suppression of symmetric and antisymmetric flutter with a 24-percent in- 
crease in attainable dynamic pressure. Posttest analyses are shown which 
illustrate the problems encountered with the other laws. 


Introduction 

Aircraft designs that emphasize the reduction of 
structural weight to maximize efficiency and agility 
increase the likelihood that active flutter suppres- 
sion will be needed to remove structural dynamic 
instabilities. In such cases, active flutter suppres- 
sion can potentially enable achievement of enhanced 
performance with lower weight. Developing methods 
to suppress flutter and reduce structural loads has 
been an objective of the active flexible wing (AFW) 
program. 

Active controls, with flutter suppression as a spe- 
cific example, are recognized to provide maximum 
performance benefits when their impact is considered 
early in the aircraft design process. Providing a com- 
plete summary of the current status of active controls 
research is beyond the scope of this paper; however, 
selected references are cited which are representative 
of published work in the areas of flutter suppression 
and gust load alleviation (refs, 1 to 29). A refer- 
ence describing the historical development of flutter 
research is also cited (ref. 30). 

This paper describes the design, test, and eval- 
uation of three flutter suppression control laws. 
The designs were part of a joint effort by Langley 
Research Center and Rockwell International Cor- 
poration to validate analysis and synthesis metho- 
dologies through the development of digital 
multi-input /multi-output control laws for an 


aeroelastic wind-tunnel model (refs. 31 and 32). The 
test vehicle used in this effort is the Rockwell AFW 
wind-tunnel model (ref. 33), which was modified from 
its initial configuration through the use of wing tip 
stores containing destabilizing mass ballast. The test 
results described in this paper refer to a model entry 
in the Langley Transonic Dynamics Tunnel (TDT) 
in November 1989. A subsequent test was per- 
formed in March 1991 during which four flutter sup- 
pression control laws were successfully tested during 
steady flight and while performing aggressive rolling 
maneuvers (refs. 34 to 37). 

The wind-tunnel model, test operating con- 
straints, tunnel turbulence model, uncontrolled flut- 
ter characteristics, and controller design constraints 
are described to provide an understanding of the con- 
troller design problem. The design objective is to 
control flutter to the tunnel limit without saturating 
control power capabilities. 

Three controllers have been designed, all of which 
use multiple control surfaces and/or sensors. The 
first controller uses traditional tools (pole/zero loci 
and Nyquist diagrams, ref. 38); the second uses 
accelerometer output blending and control com- 
mand distribution to obtain an estimate of the flut- 
ter mode rate for feedback (refs. 11, 19, and 25). 
The third controller uses linear quadratic Gaussian 
(LQG)/loop transfer recovery (ref. 39) plus order re- 
duction (ref. 40). The design approach for each is 
described, and predicted controller performance is 



shown. Test 

results are also discussed which show 

Ki,K 2 ,K 3 

control law transfer matrices (see 

that only the traditional pole/zero loci design was 
sufficiently robust to modeling errors to suppress flut- 


eqs. (1), (4), and (9)) 

ter during the test. Posttest analyses are presented 

k 

gain in desired transfer function at 

which explain the problems that were encountered 


single strand point (see eq. (6)) 

with the other laws. 

Symbols 

k \ , k‘) 
k~ 

scalar feedback gains in control laws 
(see cqs. (1) and (4)) 

smallest magnitude intersection of 

a 

parameter in washout filter element 
in control laws 


loop transfer function with negative 
real axis to left of the —1 point 
(see fig. 11); the larger fc_ is the 

a i , (12 

filter parameters in modal rate 
feedback control law (see fig. 9) 


more tolerant the system is to gain 
decrease 

j ) C dj i ) 

parameters in rational transfer 
function representation of j th 
actuator 

k+ 

largest magnitude intersection of 
loop transfer function with negative 
real axis to right of the —1 point 
(see fig. 11); the smaller k+ is the 
more tolerant the system is to gain 
increase 

b s 

wing semispan of wind-tunnel 
model 


Bln 
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control law feedthrough term (see 
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control law system matrix (sec 
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control law matrix relating outputs 
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Laplace variable 

time between samples 

matrix defining dynamic portion 
of modal rate feedback control law 
(sec fig. 9) 

analytically predicted loop transfer 
function at single-strand point in 
modal rate feedback control law 
(see eq. (5)) 

rational transfer function approxi- 
mation of j th actuator 

pseudo-control command (see fig. 9) 

output of digital control law at fcth 
time step (see eq. (2)) 

fluid velocity 

weight associated with violation of 
mt.h actuator rate constraint (see 
eq. (7)) 

weight associated with cost J (see 
eq. (7)) 

gust velocity 

general and reference streamwise 
coordinate 

discrete controller state at A:th time 
step 

continuous controller state for LQG 
law (see eq. (9)) 

spanwise coordinate 

input to digital control law at A*th 
time step (see eq. (2)) 

vector of accelerometer outputs 

achieved and desired flutter modal 
rate (sec figs. 9 and 10) 

incremental angle of attack due to 
turbulence 


6 vector of control surface deflections 
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damping ratio 

V 

white-noise input into Dryden Hit c^r 

a 

singular value 

a US ’ a <h\ 

rms symmetric and antisymmetric 
turbulence velocities, respectively 

r 

time delay 


Dryden turbulence power spectral 

density 

UJ 

frequency, rad/sec 

UJy 

break frequency in Dryden turbu- 

lence representation 

Subscripts 


A 

antisymmetric 

c 

commanded 

D 

denominator 

d 

digital representation of controller 

f 

flutter 

n 

natural frequency 

N 

numerator 

S 

symmetric 

Notation: 

Bold s^ 

/mbols refer to matrix or vector quantitite 
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is a minimum of ( ) 
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bar over symbol indicates that it is 
a maximum of ( ) 
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dot over symbol indicates time 
derivative of ( ) 

Abbreviations: 

A 

analog 

AAF 

antialiasing filter 

AFS 

active flutter suppression 

AFW 

active flexible wing 

BRF 

band-rejection filter (see eq. (8)) 

conj 

conjugate 

CPE 

controller performance evaluation 

D 

digital 
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dof degree of freedom 

DNOT notch filter to be implemented 

digitally (see “Modal Rate Feedback 
Control Law Definition” section of 
appendix) 

INTI, INT2 stable integrators (see “Modal Rate 
Feedback Control Law Definition” 
section of appendix) 

LEI leading edge inboard 

LEO leading edge outboard 

LQG linear quadratic Gaussian 

max maximum 

MIMO multi-input/multi-output 

min minimum 

MISO multi-input/single output 

PSD power spectral density 

rms root mean square 

SISO single input/single output 

TDT Transonic Dynamics Tunnel 

TEI trailing edge inboard 

TEO trailing edge outboard 

TIP wing tip 

WOF washout filter (see “Modal Rate 

Feedback Control Law Definition” 
section of appendix) 

AFW Wind-Tunnel Model and Test 
Conditions 

Wind-Tunnel Model 

The AFW wind-tunnel model is a full-span, sting- 
mounted, wind-tunnel model that can roll about the 
sting axis (fig. 1). For the flutter suppression testing 
described herein (the TDT entry in November 1989), 
the AFW wind-tunnel model was prevented from 
rolling by applying a locking pin. The model has a 
six-degree-of-freedom force and moment balance on 
the load path to the sting and an actuator that can be 
used to adjust the model angle of attack. Four pairs 
of control surfaces exist with hinge lines near the 
one-quarter- or three-quarter-chord locations. The 
actuators for the control surfaces and for the angle- 
of-attack adjustment are powered by an onboard 
hydraulic system. The model fuselage is more rigid 
than the wings. However, the sting undergoes small 


vertical and lateral translations as well as angular 
twisting about the sting axis. 



L-89-12,445 


Figure 1. Active flexible wing (AFW) model in wind-tunnel 

test section. 

Each of the three control law designs used a subset 
of four pairs of accelerometers and four pairs of con- 
trol surfaces (fig. 2). Strain gauges on the wing were 
also available, but they were not used for flutter sup- 
pression during the 1989 tests. A digital computer 
implemented the controller that processed the sig- 
nals from the accelerometers to generate commands 
for the control surface actuators; this process actively 
suppressed flutter. Signal transmission between the 
digital controller and the analog sensors and actua- 
tors required low-pass analog filters (to reduce alias- 
ing) as well as analog-to-digital and digital-to-analog 
conversions. 



Figure 2. Sketch of AFW wind-tunnel model. 


The original configuration of the AFW wind- 
tunnel model was previously tested using multiple 
control surfaces (ref. 41) to study rapid rolling ma- 
neuvers for a model with a soft, flexible wing. This 
configuration did not flutter within the operating 
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range of the TDT. For the study described here, 
tip stores containing inass ballast were added to the 
AFW. The ballast lowered the frequency of the first 
torsion mode (bringing it closer to the frequency of 
the first bending mode) and, thereby, reduced the 
dynamic pressure at which flutter occurs to within 
the operational range of the wind tunnel. 

The tip ballast stores, which are normally coupled 
in torsion with the wing via a hydraulic brake, can 
be decoupled by releasing the brake, thereby leaving 
the stores restrained in torsion by only a soft spring. 
Upon brake release, the decoupled configuration has 
flutter-free characteristics to a much higher dynamic 
pressure than those shown in figure 3 for the coupled 
case. Thus, the tip ballast stores utilize the decoupler 
pylon concept (ref. 42) to also provide a flutter- 
stopper capability. A more complete discussion of 
the tip ballast stores can be found in reference 43. 

Wind-Tunnel Conditions 

The Langley TDT, which is specially configured 
for testing aeroelastie models (ref. 44), is a sealable 
wind tunnel in which Mach number and dynamic 
pressure can be varied independently bv changing 
motor speed while simultaneously changing stagna- 
tion pressure in the tunnel through the use of pumps. 
For the 1989 test, air was the test medium, and the 
tunnel was operated in an unsealed condition at at- 
mospheric pressure. The maximum dynamic pres- 
sure achievable in this mode was variable because 
it was dependent upon atmospheric conditions. At 
standard atmospheric pressure, a maximum dynamic 
pressure? of approximately 325 psf could be generated 
at a Mach number of 0.5 (fig. 3). Higher Mach num- 
bers required air to be pumped from the tunnel to 
reduce the stagnation pressure. 



Figure 3. Wind-tunnel test conditions with air as test 
medium and predicted flutter characteristics for coupled 
configuration. 


pressure and allowed rapid progression between test 
points. Each closed-loop test, run was made at 
atmospheric pressure and was planned to reach 
a Mach number of 0.5 at the maximum achiev- 
able dynamic pressure (fig. 3). On the test path, 
both Mach number and dynamic pressure were si- 
multaneously varied by increasing the fan motor 
speed. The Mach number variation was accept- 
able because Mach number effects in the low sub- 
sonic region were small. A rapid progression be- 
tween test points was made possible by operating in 
this mode. Running at a fixed Mach number, which 
would involve pumping air out or bleeding air into 
a sealed tunnel, would have been less efficient, for 
obtaining closed-loop data. 

Wind-tunnel turbulence had a significant impact 
on the AFW wind-tunnel model response. Conse- 
quently, critical loads were monitored during the test, 
and test runs were terminated when a load exceeded 
a preset maximum amplitude. Prior to the test and 
based upon earlier tunnel entries, the turbulence was 
estimated to have a root-mean-square (mis) velocity 
of 1 ft/sec, a maximum intensity at 10 Hz. and an 
apportionment that was 85 percent symmetric and 
15 percent antisymmetric. 

Mathematical Modeling 

Linear aeroelastie descriptions for the symmetric 
and antisymmetric boundary conditions were gen- 
erated using the Interaction of Structures, Aero- 
dynamics, and Controls (ISAC) system of programs 
(ref. 45). The doublet lattice aerodynamic the- 
ory was used (ref. 46). In addition, those models 
were combined with empirical data to form a whole- 
aircraft model for a batch simulation (ref. 47) which 
considered quantization, computational delays, actu- 
ator position and rate limits, and asymmetries caused 
by differences between individual actuators. 

Actuators 

Frequency responses for the eight individual ac- 
tuators were measured with no air flow and with 
the wing elastic motion restrained. In the frequency 
range of interest , third-order transfer functions, with 
parameters optimized in a least -squares sense, pro- 
duced good fits with the measured frequency re- 
sponse data. In general, right arid left members of an 
actuator pair required different parameters to achieve 
a good fit, and, therefore, they were modeled in this 
manner. All the actuator transfer functions had the 
following form: 


Sj _ k 'h%^ dj 


j~ = 7 w, - 7 rv = T dj {> s) 

Cj (* + % ) ( s + 2 ^j s + T ) 


A test plan was developed which maximized the 
demonstrable increase in closed-loop flutter dynamic 
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where k ( [. is the steady-state gain, a c [. is the first- 
order pole location. Q is the damping ratio of the 
complex pair, and u ;</. is the natural frequency of 
the complex pair. The physical origin of these pa- 
rameters is explained in reference 47. For linear 
analyses, averages of each parameter from members 
of a pair were employed to represent each member 
to retain decoupling between symmetric and anti- 
symmetric degrees of freedom. However, for the 
batch simulation, the distinct actuator identities 
were retained and were sources of coupling between 
the two symmetries. 

Turbulence 

A basic assumption made in computing symmet- 
ric aerodynamic forces caused by turbulence was that 
randomly generated turbulence traveled, unchanged, 
downstream past the model at the fluid velocity V . 
Thus, what was encountered at a reference stream- 
wise location (x{),y) at time t would be encountered 
at a point (x, y) at a time t + (x — X{))/V , where V 
was the fluid velocity. Another assumption was that 
the symmetric turbulence down wash field seen at 
streamwise coordinate .X'o resulted in angle-of-attack 
perturbation elements of the form 

<*gs ( x 0< 2A w o, t) = exp (iu\)t) 

which had the Fourier transform 

ttf/.s- (XO. U' “0 - (w - wo) 

The same symbols (e.g., (Xy s ) arc used in this paper 
to represent a function and its transform. The 
expression for the Fourier transform at an arbitrary 
location, relative to that at the reference streamwise 
coordinate, was 

n gs (a-, y, u>) = exp ^t° ) 6(uj- w 0 ) 

The contribution, which was due to an oscillatory 
component of frequency u>o, to the antisymmetric 
turbulence was approximately proportional to the 
spanwise location, relative to the model centerline, 
with no streamwise variation. Thus, the Fourier 
transform for this case was 


The final assumptions had to do with the spectral 
content of the turbulence. A Dryden atmospheric 
turbulence model was used (ref. 48). The turbulence 
power spectral density was 


The following transfer function, used in the develop- 
ment of the state equations discussed subsequently, 
will produce an output with the power spectral den- 
sity just shown, when receiving white- noise input rj 
with the power spectral density l/n: 




1 + 3 (s 


nu>a 


1 + te) 


(0 < u < oc) 


U), 


N + 73^ 


-f (- s ) = <r w v/3 ~y 0 2 

f) S* + 2ijJgS T LJy 

which can be verified by using the following equality: 

M = - — (*w) ^ (-iuj) 

J 7T 7] X) 


A break frequency ujg/(2n) of 17.23 Hz was used to 
approximate the expected wind-tunnel turbulence. 
Resonance peaks at 10 Hz were observed in tunnel 
data from prior entries, and a range from 10 Hz 
to 12 Hz was the predicted flutter frequency. A break 
frequency of 17.23 Hz produces a peak magnitude in 
Wg/rj at 10 Hz. 

No accurate representation of the wind-tunnel 
turbulence levels was available. Furthermore, the ef- 
fect on the plant of any given turbulence intensity 
would be highly dependent on the configuration and 
the gust mode shapes used to characterize the tur- 
bulence. Based on prior wind-tunnel entries, the rms 
turbulence velocity magnitude was estimated to be 
approximately 1 ft/sec at the wing tips at a Mach 
number of 0.5 and a dynamic pressure of 300 psf. 
This magnitude was conservatively assumed to apply 
at all test conditions. Eighty-five percent of the tur- 
bulence was allocated to the symmetric component, 
and 15 percent was allocated to the antisymmetric 
component at the wing tips. 


y 

a <h\ ( x ’ y-,^) = ^-6 {u - u/’o) 

The parameter b s was the wing semispan. The coeffi- 
cients of the Dirac delta functions in these downwash 
representations were used to compute generalized 
aerodynamic forces caused by turbulence. 


In the batch simulation, a single Gaussian- 
distributed random number with unit standard de- 
viation was generated at each integration time step 
(0.5 msec) and scaled to provide a digital approxi- 
mation of continuous white noise with a power spec- 
tral density of 1 / 7r (ref. 47). That single number 
was then multiplied (by either 0.85 or 0.15) and used 
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as an input to the separate transfer functions rep- 
resenting symmetric and antisymmetric turbulence. 
Given the gust mode shapes just discussed, this pro- 
cess resulted in a constructive correlation on the 
right wing (at least at the X{) location), a destructive 
correlation on the left wing, and a linear variation 
accross the span. Using two uncorrelated random 
numbers would have been preferable with the rms 
gust velocity components of a (JS = 0.985 ft/sec and 


(T g A = 0.174 ft / sec to achieve an 85/15 distribution 
whose vector sum was 1 ft./sec. The representation 
employed was sufficient to evaluate the potential for 
saturating the actuator rate capability, because the 
larger of the values for the left and right control sur- 
face activity was used. In the results, the predicted 
actuator rate rms responses will be seen to bo sig- 
nificantly higher than those in the test actuator rate 
rms responses. 


Equations of Motion 

In vacuo vibration characteristics were' calculated based upon a structural model developed bv Rockwell 
International Corporation. The structural model development benefited from experimental data obtained for 
the model prior to the addition of the tip ballast stores. The vibration tests were 1 also made on the current 
configuration, and the results from these tests were used to adjust the predicted modal natural frequencies and, 
together with aeroelastic analyses, to determine which elastic modes to retain in the model. Eight symmetric 
and seven antisymmetric elastic modes were retained. Structural damping was modeled as being viscous with 
a damping ratio of 0.015 assumed for each mode. 

Equations of motion were developed both in a frequency domain form, which made direct use of the tabular 
unsteady aerodynamic forces, and in a finite dimensional state-space form, which employed rational function 
approximations to the aerodynamic forces. The frequency domain form of the equations was 


nij (iu>) 2 + 2Qu> nj (iu>) +uf, . - </Q‘ 


M/, (iw) 2 

I 



( l Qf} ^ / (*g 

o T h bj k \ \6r 


with accelerometer outputs of the form 


zj = (iw) 2 n(xj,yj)Z{u)) 

A display of the dependence of the unsteady aerodynamic forces upon frequency and Mach number was 
suppressed for brevity and clarity. All the mathematical models were generated with the aerodynamic force 
coefficient data corresponding to a Mach number of 0.5, regardless of velocity, because Mach number effects in 
the low subsonic region of the test were small. This approximation became more precise as the tunnel operating 
limit was approached. The row vector Tl(xj,yj) relates the jth output to unit displacements of the generalized 
coordinates £. The Kronecker delta is zero unless the indices match, in which case it is unity. This form 
of the equations allowed the numerical computation of a frequency response function for any output/input 
pair. Consequently, for a stable system, power spectral densities and rms values for any output could also be 
computed. The modal rate feedback design approach used design models with the frequency domain form. 

Rational function approximations were made to the unsteady aerodynamic forces to obtain finite dimensional 
state-space models. Reference 49 presents this approach. The same single lag factor per mode was judged 
to sufficiently represent the unsteady aerodynamic effects, while keeping the number of states manageable. 
The resulting state-space models with three pairs of actuators retained were 35th order symmetrically and 
32nd order antisymmetrically. State-space plant representations were employed in the traditional pole/zero 
loci and modified LQG design approaches. Reference 47 presents the structure of the state-space models and 
the additions necessary to generate the batch simulation. 

Controller Design Considerations in Active Flutter Suppression 

Flutter Character 

For the 1989 test described here, the model was fixed in roll. This condition, together with the nearly rigid 
fuselage, caused the characteristics of the AFW wind-tunnel model symmetric and antisymmetric flutter to 
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be similar. The predicted flutter dynamic pressures at a Mach number of 0.5 (see fig. 3) were 248 psf for the 
symmetrical and 252 psf for the antisymmetrical cases. Therefore, the simultaneous operation of a control law 
for each symmetry was required. The flutter frequency in each case was predicted to be approximately 11.5 Hz. 
The flutter characteristics were relatively explosive; at a dynamic pressure of 325 psf, the predicted time for 
the flutter mode to double in amplitude was approximately 1/10 sec. 

Control Law Implementation Considerations 

Each of the three active flutter suppression (AFS) control laws was designed in the continuous domain. 
The assumption was made that no coupling existed between the symmetric and antisymmetric response for 
the AFW wind-tunnel model. Figure 4 illustrates how the symmetric and antisymmetric forms of the control 
laws were implemented simultaneously by the digital controller (ref. 50). For each pair of accelerometers, the 
symmetric signal was determined as the average of the right and left signals, and the antisymmetric signal was 
determined as one-half of the difference between the right and left signals. Similarly, the right and left control 
surface commands were determined as the sum and difference of the symmetric and antisymmetric commands 
for each pair of control surfaces. 

The trailing-edge outboard (TEO) control surfaces (see fig. 2) were the most effective in controlling flutter, 
although the actuator hinge moment available for these surfaces was only one-half as large as the others because 
each TEO surface was driven by one actuator rather than two; the use of a single actuator was due to the 
limited space available in the outboard portion of the wing. The leading-edge-outboard (LEO) surfaces have 
unfavorable aerodynamic loading that does not restore the surfaces to a neutral position if the actuators become 
overloaded. The trailing-edge-inboard (TE1) surfaces have favorable aerodynamic loading, but they are not as 
effective as the TEO surfaces in controlling flutter. The leading-edge- inboard (LEI) control surfaces had higher 
inertias and less aerodynamic effect upon the flutter mechanism and were, therefore, not suitable for FSS con- 
trol. Each of the wing accelerometer pairs was located near the hinge line of one of the control surface pairs, with 



Figure 4. Implementation of symmetric and antisymmetric AFS control laws. 
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the exception of the wing tip (TIP) accelerometers 
that were located approximately midchord near the 
wing tips. The TIP accelerometers responded to the 
flutter mode and, at the same time, were relatively 
unresponsive to the higher frequency modes when 
compared with the other accelerometers. 

The digital implementation of the control laws 
had certain implications for the control law designer. 
The sample rate was 200 Hz. A low-pass analog 
antialiasing filter was required for each channel be- 
ing digitized to attenuate signal strength above the 
Nyquist frequency (100 Hz in this case), so that 
higher frequency harmonic signals would not corrupt 
the lower frequency signals caused by the periodic 
sampling. Two viable choices of analog antialiasing 
filters were provided. The choices were a first-order 
filter with a break frequency of 25 Hz and a fourth- 
order Butterworth filter with a break frequency of 
100 Hz. Each of these introduced approximately 
the same lag in the flutter frequency range (near 
10 Hz). Only the first-order filter was used during 
the test, principally because it provided some atten- 
uation of outputs that were due to structural modes 
in the 30- to 40-Hz range; these outputs were out- 
side the desired control bandwidth but within the 
actuator bandwidth. A first-order filter with a break 
frequency of 100 Hz and a fourth-order Butterworth 
filter with a break frequency of 25 Hz were also avail- 
able; these filters were not viable candidates because 
of the lack of antialiasing protection from the former 
and the excess phase lag from the latter. 

The signal amplitude was quantized because 
the analog-to-digital and digital- to-analog convert- 
ers had 12 bits of resolution. Using an accelerome- 
ter signal as an example and assuming that the ac- 
celerometer signals had a range of 40 g (±20$) led to 
a quantization of approximately 0.01$. This quanti- 
zation level is of little concern unless the controller 
has an extremely high gain such as would occur at 
low frequency for a pure integrator. 

An effective average 0.5 time step delay was in- 
troduced by the sampling because, after a signal was 
sampled at the beginning of a time step, no addi- 
tional information was passed until the beginning of 
the next time step. The control laws were imple- 
mented in the digital computer so that control com- 
mands generated based on sensor inputs received at 
one sampling instant were held (before being sent 
out) for a time T/, such that the sum of and r c (the 
time required to compute the outputs) was 5 msec 
(one time step). Thus, approximately a 1.5 time 
step delay was present as a result of these aspects 
of the digital implementation as compared with a 
continuous implementation. 


Each of the controllers was designed in the con- 
tinuous domain and implemented digitally by using a 
Tustin transformation with no frequency prewarping 
(ref. 51). The method chosen for transformation of a 
continuous controller to a discrete form also impacts 
the digital time delay. The net time delays assumed 
in the three designs will be shown below. 

Design Objectives and Requirements 

The design objective for all three flutter suppres- 
sion control laws was to demonstrate closed-loop sta- 
bility up to the wind-tunnel limit of 325 psf dy- 
namic pressure. This objective would constitute 
a 30-percent increase in the flutter dynamic pres- 
sure relative to the lowest predicted open-loop flutter 
boundary in the subsonic region. 

For those control laws that had a single-strand 
(i.e., single-channel) point in the feedback loop, a 
predicted gain margin of ±6 dB and a predicted 
phase margin of 4=30° with respect to that point 
were required throughout the test envelope. (This 
phase requirement was a relaxation from an earlier, 
unachieved, requirement of 45°. ) If the law was 
not truly single input/single output (SISO), then 
the simple gain and phase margins were a poten- 
tially nonconservative assessment of robustness, and 
they will be shown to be a source of poor con- 
troller performance for one of the control laws. I he 
control law that utilized multi-input /multi-output, 
(MIMO) analysis was judged by potentially conser- 
vative multi- variable margins, and the stated require- 
ments were relaxed for that case to correspond to the 
predicted robustness achieved. 

Actuator rate saturation can degrade control law 
performance and lead to closed- loop instability. At 
the wind-tunnel limit of 325 psf, the open-loop time- 
to-double amplitude for the symmetric flutter mode 
was predicted to be 1/10 sec. For this instability 
level, actuator rate saturation of a pair of actuators 
for even a brief time in response to wind-tunnel tur- 
bulence could cause an unacceptably large growth 
of the flutter mode. This potential for rate satu- 
ration reinforced the need to restrict actuator rate 
requirements. The TEO control surfaces were pre- 
dicted to have a peak no-load aerodynamic rate ca- 
pability of 225 deg/sec. It was assumed that no rate 
saturation for a three-standard-deviation turbulence 
velocity magnitude was adequate for assuring that 
rate saturation was sufficiently unlikely (ref. 52). A 
one-standard-deviation actuator rate of 75 deg/sec 
resulted. At a predicted flutter frequency of approx- 
imately 11.5 Hz, this rate constraint translated to 
a maximum rms control deflection of 1° (i.e., the 
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Figure 5. Sensors and control surfaces used by each control law. 


maximum deflection that would occur if all control 
power was concentrated at the flutter frequency). 

Another requirement for evaluating candidate 
control laws prior to the wind-tunnel entry was to 
demonstrate closed-loop stability throughout the test 
envelope using a batch simulation (ref. 47). As dis- 
cussed, the simulation replicated quantization effects 
because of finite word length in the signal convert- 
ers and imposed rate and displacement limits on the 
control surface actuators. This simulation also al- 
lowed both symmetries to be run simultaneously in 
the presence of simulated turbulence excitation with 
a separate dynamic actuator model for each left and 
right control surface. 

Design and Analysis of Active Flutter 
Suppression Controllers 

Three active control laws for flutter suppres- 
sion were designed in the continuous domain, imple- 
mented digitally, and tested in the TDT. The dark- 
ened areas in figure 5 show the sensors and control 
surfaces employed in each law. The next three sub- 
sections present the design approach and pretest 
analysis of performance for each control law. The 
approaches are presented in the order of increasing 
mathematical complexity. 


Traditional Pole/Zero Loci Design (Design 

Number 1) 

Overview . This control law was designed us- 
ing traditional complex plane mappings of poles and 
zeros. The principal philosophy behind this design 
effort was to avoid getting lost in complexities that 
arc of secondary importance with respect to the flut- 
ter control problem and to reduce the problem to its 
bare essentials. One step toward accomplishing this 
was to concentrate primarily on the two structural 
modes that participate directly in the flutter and on 
the SISO zeros in the same frequency range which 
result from the choice of a particular sensor pair and 
control surface pair. Vibration modes dominated by 
sting deflections and their associated zeros were ig- 
nored, as were high-frequency modes. This procedure 
was possible because of effective pole/zero cancella- 
tions associated witli the chosen control surfaces and 
sensors and because of frequency separation between 
the flutter dynamics and high-frequency modes. 

Selection of sensors and control surfaces was a 
necessary first step in the controller design. The 
accelerometer pair at the TIP locations was chosen 
because this pair was responsive to the flutter motion 
and least responsive to high-frequency modes. The 
fourth-order 100-Hz Butterworth filter was used in 
the controller design and analysis. The TEO control 
surface pair was chosen because it was the most 
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effective in controlling flutter. The TEI control 
surface pair was added later in the antisymmetric 
control law to improve stability margins. The design 
plant, and final controller structure are shown in 
figure 6. In equation form, the controller was 

6 r = K ] (s)z vl v 


d 

1 


k\ 


fi 

s + a 


• s ^ + 2Cv^.V' s + 

+ 2(i)ujps + u/jj 


2TIP 


(1) 


where 


Sr 


^TEI, 

^TEOr 


The symmetric law was SISO. (The switch shown in 
fig. 6 was open.) All parameters in the controller 
were the same for both symmetries, and they have 
been specified numerically in the appendix. The term 
K t (s) was implemented digitally for testing in the 
TDT. 


Design plant, (s) 



Compensator dynamics 


Figure 6. Block diagram of traditional pole/zero loci control 
law. 


The steps to arrive at this controller form are 
presented subsequently. Note in figure 6 that the 
commands sent to the TEO and TEI surfaces in the 
antisymmetric case were dynamically equivalent in 
that they differed only by the constant factor ( d ). 
SISO design and analysis techniques were employed 
despite the possibility, for the antisymmetric law, 
that the resulting robustness characterization was 
nonconservative. 

Straight feedback with no dynamic compensation 
was investigated first for an SISO design using the 
TIP accelerometer pair and the TEO control surface 
pair to see whether this feedback would be sufficient 


to stabilize the system, and if not, what problems 
would be encountered when attempting to employ a 
simple solution. Consideration of the high gain that 
was required and the desire to ensure a favorable root 
locus path led to the use of a second-order inverted 
notch filter to be described later in this section. 
A final consideration was that the response of the 
system to steady-state bias errors must be acceptably 
small: this consideration led to the addit ion of a first- 
order washout filter and resulted in a third-order 
controller. 

Critical zero . For t he sketches of poles and zeros 
presented in figure 7, the horizontal axis was greatly 
exaggerated relative to t he vert ical axis to show more 
detail. All poles and zeros not. associated with the 
compensator should be considered to lie near the 
imaginary axis. The sketch in figure 7(a) shows 
the loci of poles and zeros as functions of dynamic 
pressure. The poles are those associated with the 
two strongly interacting modes for the AFW wind- 
tunnel model wit h no active compensation, and the 
zeros arise from a particular choice of sensors and 
actuators. The pair of superimposed zeros at the 
origin results from the fact that accelerometers were 
used for feedback. A critical zero closely associated 
with the higher frequency of the two interacting 
modes was found for the TEO control surface and 
the TIP sensor. As dynamic pressure increased, the 
critical zero and the pole associated with the higher 
frequency mode stayed near each other until just 
below the flutter dynamic pressure at which point 
the pole broke away to the right and crossed into the 
unstable, right half of the complex plane. 

The use of simple feedback will drive the closed- 
loop roots from the open-loop poles to the transfer 
function zeros as a function of feedback gain. How- 
ever, given uncertainties in the model of the plant, 
it is not always clear what path the roots will take 
(ref. 27). Figure 7(b) shows how the system might 
be stabilized by simple feedback; figure 7(c) shows 
a case in which no value of feedback gain exists for 
which the closed-loop system will be stable. 

Even when the desired path is followed, the loca- 
tion of the critical zero near the imaginary axis indi- 
cates that a high gain would be required to drive the 
unstable root close to the zero to stabilize the system. 
One difficulty associated with high-gain controllers is 
that the control surface rates required to control flut- 
ter while subject to continual turbulence excitation 
would be large and might saturate the rate capabil- 
ity of the actuators, thus causing a loss of flutter 
control. Another difficulty associated with high-gain 
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Pi = First structural mode 
p 2 = Second structural mode 


z c = Critical zero 


(a) Open-loop poles and zeros as function of dynamic pressure. 



(b) Stabilizing feedback gain loci. 



(c) Nonstabilizing feedback gain loci. 

Imaginary 



Pp = Filter pole 
Zp = Filter zero 


(d) Stabilizing dynamic feedback gain loci. 

Figure 7. Sketch of pole/zero loci. Arrows indicate direction of increase in either dynamic pressure or gain. 


controllers is that high-frequency modes or actuator 
roots can be driven to be unstable (ref. 53). 

Inverted notch filter . Dynamic filtering was 
required to reduce the feedback gain required for 
stabilization because of the location of the critical 
zero near the imaginary axis. The intent for this 
control law was to “soften” the effect of the critical 
zero by placing a filter pole near the critical zero and 
placing a filter zero farther to the left, as shown in 
figure 7(d). The result is similar to an inverted notch. 

The location of the critical zero changes as a func- 
tion of dynamic pressure, whereas the location of 
the open- loop filter pole is independent of dynamic 
pressure unless scheduling of controller dynamic 
parameters is used. Also, the actual locations of 


system zeros are difficult to predict analytically and 
can be difficult to measure experimentally. To avoid 
scheduling and because of the uncertainty about the 
exact location of the critical zero, the filter pole was 
placed somewhat to the left of the predicted criti- 
cal zero and had a damping ratio of approximately 
10 percent. This placement assured that the result- 
ing pole-zero interaction would cause the desired sta- 
bilizing root- locus path to be achieved, even when 
subjected to moderate plant variations and model- 
ing errors. Using frequency domain Nyquist crite- 
ria for stability margin analysis, a 50-percent damp- 
ing ratio at a natural frequency 20 percent higher 
than that for the compensator pole was chosen for 
the compensator zero. This zero choice was made to 
simultaneously maximize the gain and phase margins 
of the system over a range of dynamic pressures. 


12 


Figure 8 shows the magnitude and phase plots of 
the frequency response of the inverted notch filter. 
The filter amplifies the control surface activity in the 
frequency range predicted for flutter, which in this 
case is approximately 11.5 Hz. Because the control 
surface activity was concentrated at this frequency, 
the controller made efficient use of the available 
control power and was fairly insensitive to modeling 
errors outside the frequency range of interest. 


O) 



Frequency, Hz 

Figure 8. Frequency response for inverted notch filter in 

pole/zero loci control law. 

It is generally thought that successful designs us- 
ing pole/zero cancellation require an accurate knowl- 
edge of the plant; however, the inverted notch filter 
was evaluated using variations in the model of the 
AFW. The control law was judged to be tolerant to 
changes in the frequency of the flutter mode, and this 
tolerance was due in part to the robust placement of 
the filter pole with respect to the critical zero. The 
result was that the stabilizing character of the root 
locus did not change despite the frequency shifts. 

The numerical values of the controller parameters 
in the continuous domain arc presented in the ap- 
pendix. The controller was implemented digitally by 
using a Tustin transformation and a 200-Hz sample 
rate. The design plant employed by the pole/zero loci 
methodology did not account for the time delays that 
resulted because of the digital implementation (see 
fig. 6). A “buy-back” approach was implemented, 
however, which approximately removed the effect of 
the 1 time step delay that was a result of the way 
the control law was implemented in the digital com- 
puter. The controller implementation was such that 
the controller output, which was computed based on 
sensor inputs and states at time was not sent out 
until The buy-back procedure follows. The 


discrete state equations for the controller resulting 
from the Tustin transformation were 

x (A'+i) = f </ x a + G,m- 1 
«A = H ,/x*. + E,/y* / 

or, at t {k+l] 

U (A+ 1 ) = H f/ X (A+1) + E -/y(A+l) 

Because of the time delay in sending out the con- 
trol command, the following was implemented if the 
buy-back procedure was not employed: 

u (A+i) = H f ,x, + E <t y k . 

In the buy-back procedure, the designer replaces H ( y 
with = H,/F f y and E f y with = (H,yG f y + E,y) 
in the implemented equation and makes the approx- 
imation that yj. ~ One obtains 

u (A+l) = HW + E^y k . 

= H,/ x (t-+i) + E,/y*. 

~ H ,/ x (fr+i) + E </y(A+i) ( :i ) 

The smaller that E,yy^, +1 j and E ,/y^. are (relative to 
H ( yx^._j_^) the better the approximation. The ap- 
proximation is exact if E,y — 0. The lead int roduced 
provides an approximate buy-back of the one-step 
computational delay. 

Predicted performance . The traditional 

pole/zero loci design was predicted through linear 
analysis and batch simulation to provide closed-loop 
stability up to the limit of the wind-tunnel operating 
range, as summarized in column two of table I. The 
term r/ max refers to the maximum dynamic pressure, 
measured in pounds per square foot, for which the 
closed-loop system was predicted to be stable. 


Table 1. Predicted Performance for Traditional 
Pole/Zero Loci Design 




Margins 
at 300 psf 

mis cont rol 
activity 
at 300 psf 
(percent of 
max allowed) 

Degrees 

<7iuax • 

Gain. 

Phase. 

TEO, 

TEI. 

of freedom 

psf 

clB 

<le K 

percent 

percent 

Symmetric 

> 350 

±7 

±33 



Antisymmetric 

>350 

±7 

±38 

t «> 

*5 
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A common evaluation point of 300 psf was chosen 
for each of the control law's. The gain and phase mar- 
gins for the pole/zero design w r ere predicted through 
linear analysis to exceed the design requirements at 
300 psf. Other analyses, not showm, predicted that 
the margins were maintained throughout the test 
envelope. Positive and negative gain margins w'ere 
verified in the batch simulation at selected dynamic 
pressures by individual variation of symmetric and 
antisymmetric gains until simulation time histories 
exhibited oscillatory divergence. Phase margins were 
more difficult to verify in the batch simulation, and 
only tolerance to phase lag was investigated. This 
investigation was accomplished by incrementally re- 
ducing the break frequency of the fourth-order But- 
terworth filter from the nominal value of 100 Hz until 
the simulation time histories show r ed oscillatory di- 
vergence. The tolerance to phase lag was equated 
to the additional phase lag because of the break 
frequency reduction at the frequency of divergent 
oscillation. 


The rms control activity in table I is shown in 
terms of a percent of the maximum acceptable rms 
rate of 75 deg/sec. The predicted requirements 
are well within their limits, and the TEO surfaces 
are the dominant ones used for flutter suppression. 
These results have been generated using the batch 
simulation in the presence of simultaneous symmetric 
and antisymmetric turbulence excitation. 


Analyses were made which predicted that the 
pole/zero loci control law structure was robust to 
variations in the plant. In these analyses, an early 
version of the inverted notch filter was employed 
for which the compensator zero w T as more heavily 
damped ( (\\ T = 0.707) than that ultimately selected 
(<^y = 0.47). The robustness characteristics, to be 
described, hold for both choices; the primary reason 
for the ultimate choice of the less-damped zero was 
to simultaneously maximize the gain and phase mar- 
gins of the system over a range of dynamic pressures. 
Cases were examined with the SISO form (the swatch 
showm in fig. 6 w r as open) with both symmetries; and 
cases w^ere examined with the MISO form (the swatch 
was closed). Each version of the controller stabilized 
the plant to a dynamic pressure of 325 psf or more un- 
der the following simulated conditions: a Mach num- 
ber of 0.5 in air, a Mach number of 0.8 in a heavy gas 
test medium, and a Mach number of 0.9 in a heavy 
gas test medium. Because the controller stabilized 
these varieties of plants, there was confidence that 
wind-tunnel testing of the pole/zcro loci controller 
would be successful. 


Modal Rate Feedback Design (Design 

Number 2) 

Overview. The design philosophy for the modal 
rate feedback control law' was to use linear combina- 
tions of multiple accelerometer signals and dynamic 
compensation to synthesize the flutter mode rate for 
feedback to multiple control surfaces (refs. 11, 19, 
and 25). The control structure used is shown in fig- 
ure 9. A 1.5 time step delay was included in the de- 
sign plant to account for the effects of a digital imple- 
mentation of a continuous controller. Consequently, 
the continuous controller had a 1.5 time step lead 
to counterbalance the digital implementation delay. 
The first-order 25-Hz antialiasing filters were selected 
for analysis and design. Numerical values determined 
for the controller parameters are presented in the ap- 
pendix. Multiple sensors were used to identify the 
activity of the flutter mode not only by frequency 
but also by the geometry of its characteristic mode 
shape. Multiple control surfaces were used to control 
the flutter mode without affecting other modes. Iso- 
lation of the flutter mode was determined by the fil- 
tering included in the design plant, the compensator 
dynamics, and the extent to which the blending and 
distributing rejected feedback interaction with other 
modes. 

Control structure definition . The blending 
coefficients used for the accelerometer pairs, the dis- 
tribution coefficients used for the control surface 
pairs, and the overall system gain were the design 
variables used in an optimization procedure. The 
ordering of the four pairs of sensors used was 


^LEO 


^TEI 
^ TEO 


L ztip J 


and the ordering of the two pairs of control surfaces 
used was 


S c = 


^TElc 
. 6 TEO,- . 


In equation form, the control law was 


6 r = K 2 (S)Z 

= Dis*2T 2 (s) B ln z (4) 

The controller K 2 (s) was implemented digitally for 
testing in the TDT. 

The controller dynamics T‘ 2 ($) were chosen by 
the control law^ designer. In the flutter frequency 
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Design plant P 2 (s) 



range, rate and position components are included in 
the single-strand signal by passing signals D\ and B 2 
through one and two integrations, respectively. Low- 
pass filter elements having break frequencies that 
were low relative to the flutter frequency, referred to 
here as stable integrators, were used in lieu of pure in- 
tegrators to avoid potential high-gain difficulties with 
sensor bias and quantization effects. The second sta- 
ble integrator (position path) was included to allow 
more freedom to change the phasing of individual 
sensor channels. A first-order “washout” filter with a 
zero at the origin and a pole at a frequency below the 
flutter frequency was also used to reduce response to 
bias errors and low-frequency disturbances. Finally, 
a second-order notch filter, which was ultimately im- 
plemented digitally, was used to adjust the phasing 
of the control action at the flutter frequency and to 
reduce the response to a nonflutter mode. The notch 
for the symmetric case was for a 5.7-Hz sting mode, 
and the notch for the antisymmetric case was for an 
18.3-Hz structural mode. 

Coefficients and parameter values are shown in 
the appendix; the process by which the design vari- 
ables were determined will be described. When 
examining the blending coefficients, note that the 
integrators attenuate a signal proportional to the 
frequency of the signal. Therefore, to more readily 


compare the contributions of individual sensors list'd 
in feedback, each row of the blending coefficient- ma- 
trix was separated into a multiplicative factor and 
that row normalized by the factor. The normalizing 
factor chosen for the first row (one integration) was 
70 rad/sec (11.14 Hz), which is a number close to the 
predicted flutter frequency. The normalizing factor 
chosen for the second row (two integrations) was the 
square of that for the first row. 

The method compares the predicted and a refer- 
ence desired response (hiring the opt imization. Thus, 
for each symmetry, individual frequency responses 
were required for each (accelerometer pair/ actuator 
pair). These responses can be obtained either 
through pretest modeling or through experiment as 
was done successfully in reference 54. Analytical fre- 
quency responses were generated using the IS AC sys- 
tem of programs (ref. 45). The responses were com- 
puted using the nonstate-space form of the equations 
of motion, thus removing the need to make ratio- 
nal function approximations of the unsteady aerody- 
namic forces. Experimentally derived frequency re- 
sponses, which were also used during the wind-tunnel 
entry to improve the control law, will be discussed 
subsequently. 

Flow diagram. The signals designated in fig- 
ure 9 as 6 C and z represent, respectively, commanded 
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control surface deflections and measured local ac- 
celerations for a given symmetry. The frequency 
response for each acceleration signal resulting from 
each control surface command (either symmetric or 
antisymmetric) for continuous analytical models of 
the AFW wind-tunnel model was precomputed for 
the frequency range from 2 Hz to 64 Hz; these compu- 
tations were repeated for several dynamic pressures 
and the frequency responses were retained for further 
analysis. For the purpose of control law design 
(performed in the continuous domain), frequency re- 
sponses representing the effects of the time delays (a 
1.5 T delay was assumed) and of a candidate set of 
analog antialiasing and notch filters were also com- 
puted and combined with the frequency responses 
that represented the AFW wind-tunnel model. This 
design plant model was used in the development of 
the control law. Frequency responses representing 
the two paths through the controller dynamics were 
also precomputed and stored in combination with 
frequency responses for the design plant. This pro- 
cess reduced the amount of time required for each 
iteration of the optimization. 

Optimization strategy. Although the normal- 
ized blending and distribution matrices contained 
eight and two coefficients, respectively, the magni- 
tude of the largest coefficient of each matrix was 
factored into a system gain. Thus, these matrices 
had seven and one degrees of freedom, respectively, 
which, together with the system gain, constituted 
the nine available degrees of freedom (identified in 
fig. 9) for use as design variables in the optimization 
procedure. 

The output from the dynamic compensator was 
a single-strand point for the feedback path. The 
optimizer was used to drive the composite frequency 
response at that point to match a simple, desired 
frequency response /?.(.s) (see eq. (6)), which was 
proportional to the idealized modal rate of the flutter 
mode, generated from a continuous state-space model 
of the AFW. Thus, the desired frequency response 
was proportional to that of the rate of a simple 
oscillator with damping ratio £/• At the design 
point chosen, 325 psf, this oscillator was unstable. 
Figure 10 illustrates the concept of feedback of modal 
rate to add damping. The value of the gain ( k ) 
was chosen to be |4(/| where the damping ratio was 
that predicted for the design point of 325 psf. For 
an unstable plant, this choice of gain corresponds 
to what would be found for the minimum energy 
stabilizing controller. 

Figure 11 shows Nyquist (or polar) plots of the 
predicted T^(s) and desired R(s) loop-transfcr- 


function frequency responses, as defined at the single- 
strand point (see fig. 9 and eqs. (5) and (6)): 

T SS (.9) = T 2 (s) B ln P 2 (s) D is fc 2 (5) 


and 


R(s) = 


kujf 


s 2 + 2 ( fUJjs + 




(6) 
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Figure 10. Block diagram of rate feedback. 



GM = r- 20l °9lO k - 1 
\-20log 10 k + / 

Figure 11. Predicted and desired single- strand loop transfer 
function; antisymmetric, 325-psf condition. 


The frequencies used in the analysis range from 
2 Hz to 64 Hz. A full Nyquist plot would span 
frequencies from minus infinity to plus infinity and 
be a symmetric function of frequency (with respect 
to the real axis) so that the depiction of the neg- 
ative frequency portion would be redundant. The 
Nyquist stability criterion requires that,, for each un- 
stable pole of the open-loop system, the Nyquist plot 
must form one counterclockwise encirclement of the 
— 1 point in order for the closed-loop system to be sta- 
ble. For oscillatory instabilities, the unstable poles 
occur in complex conjugate pairs, thus requiring two 
encirclements per pair. However, one of the encir- 
clements would occur for the frequency range from 
minus infinity to zero, which is not shown. Fig- 
ure 11 represents a condition considerably above the 
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predicted flutter dynamic pressure, and the counter- 
clockwise encirclement of the —1 point indicates 
that the unstable flutter mode would be stabilized 
through feedback. 

For actively stabilized flutter, the encirclements 
will occur in the vicinity of the flutter frequency. 
Gain and phase margins, with respect to errors at 
the single-strand point, can be read directly from a 
Nyquist plot (as indicated in fig. 11) as the amount 
of shift that can be tolerated while still encircling 
the —1 point. Errors that contribute to excess phase 
lag at the flutter frequency wdll shift the positive 
frequency plot clockwise (and the negative frequency 
plot counterclockwise) until closed-loop instability 
is encountered at a frequency slightly above the 
flutter frequency. Similarly, errors that contribute 
to excess phase lead result in closed-loop instability 
at a frequency slightly below the flutter frequency. 

The response of modes other than the flutter 
mode will be evident as additional “lobes 1 ' on the 
Nyquist plot. To the extent that the sensor and 
control surface blending can isolate the flutter mode, 
these extra lobes will be small. If these lobes are not 
small, they could result in clockwise encirclements 
of the —1 point, thus indicating that an open-loop- 
stable mode would be driven unstable, through feed- 
back, at a frequency other than the flutter frequency. 

The cost function for the optimization contains 
the sum of the squares of the difference between the 
predicted and the desired response, multiplied by 
frequency dependent weights, together with penalty 
contributions for rms actuator rate violations; this 
function is shown in equation (7): 

J = [7’ss ( ) - /?bW)] t-onj [T ss - K{ r^» )} 

n= 1 

+ [V *rn { miiX [ () ’ ( < W)ni,s " (W)r„« ] } ( 7 ) 

m=l 

where R(iu> n ) is the desired modal rate feedback fre- 
quency response at the nth frequency; T$$(iuj n ) is 
the analytically achieved single-strand frequency re- 
sponse at the nth frequency; W^ n is the weight defin- 
ing the cost of error between desired and achieved 
single-strand frequency responses at the nth fre- 
quency; (^m r ) rillH i s the maximum allowable com- 
manded rms rate for the mth actuator; (<5 mr ) nns is 
the commanded closed-loop rms control rate in re- 
sponse to turbulence for the mth actuator; and 
is the weight penalty for violating the mth constraint. 
To the extent that the actual response matches the 
desired response, the system, as observed at the 


single-strand point, will behave as though it had rate 
feedback for adding damping to the unstable flutter 
mode. 

A Davidon-Flctcher- Powell optimization routine 
(ref. 55) was used to find the system gain and blend- 
ing and distribution matrix coefficients for which 
the cost function was minimized at a particular dy- 
namic pressure. Each resulting design was evalu- 
ated at other dynamic pressures. Predicted perfor- 
mance was satisfactory throughout the wind-tunnel 
test envelope. 

Several schemes for selecting the frequency- 
dependent weights were considered. Initially, the 
weights were chosen to more heavily penalize the er- 
rors in a discrete band about the flutter frequency 
than those outside the discrete band. Alternately, 
the weights were chosen which were proportional to 
the magnitude of the desired response, thus again 
most heavily penalizing errors at the center frequency 
of the flutter mode. Uniform weighting across a lin- 
ear distribution of all available frequencies, from 2 Hz 
to 64 Hz, was also used. The weighting used for the 
final analytically and experimentally derived control 
laws was uniform on a logarithmic frequency scale, 
which penalized low-frequency errors more heavily 
than high-frequency errors compared with the linear 
frequency distribution. The uniform weighting (on 
a logarithmic scale) scheme exploits sensor blending 
and control distribution in the low-frequency range 
where the model is expected to be more accurate'. 
The scheme requires reliance, instead, upon dynamic 
filtering to attenuate' model response at high frequen- 
cies where the analytical model is less well-known. 
In practice, achievement of the desired Nyquist plot 
for the single-strand point was not particularly sen- 
sitive' to the weighting scheme chosen. Throughout 
the design process, the init ial part of an optimization 
run would capture the bulk of the achievable' Nyquist 
plot shaping at the single-strand point, with only 
limited and very slow improvement in the cost func- 
tion (and Nyquist plot) as the optimization was al- 
lowod to continue. However, none e)f the convergence? 
criteria were? achieved, and the values of the blend- 
ing and distribution coefficients changed over time, 
thus suggesting that multiple combinations of coef- 
ficients existed which would generate approximately 
the same single-strand Nyquist plot and cost func- 
tion value. Additional discussion e>f problems with 
this approach will be presented in the section entitled 
“Post test Analysis. 11 

At all times during the optimization process, 
the distribution coefficient for the TEO surface 1 was 
larger in magnitude than the distribution coefficient 
for the TEI surface. At some times, the coefficient 
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for the TEI surface had a sign opposite to the sign 
of the coefficient for the TEO surface. Opposing 
signs are reasonable for a condition in which a node 
line for the flutter mode is located between the two 
control surfaces. However, the closed-loop rms con- 
trol rate because of turbulence for the TEO surface 
was found to be higher for a suboptimal solution 
with opposing signs than for a suboptimal solution 
with like signs for the two surfaces. The opposition 
of signs occurred after a lengthy optimization cycle 
that only slightly improved the cost function and the 
single-strand Nyquist plot. (The rms control rate 
constraints were not active at this time.) Because 
of concern that the two surfaces might be “working 
against each other,’ 1 the sign of the TEI surface was 
forced to match the sign of the TEO surface, and the 
magnitude of the coefficient for the TEI surface was 
set to a value that reduced the rms activity of the 
TEO surface. This procedure removed one degree of 
freedom in the optimization. 

Control surface rates . Figure 12 shows the 
predicted power spectral density (PSD) plots at a 
dynamic pressure above the open-loop flutter point 
for symmetric closed-loop rates for the TEO and 
TEI control surfaces caused by the modeled tur- 
bulence. (A 1 -ft /sec rms, symmetric gust velocity 
was assumed.) These plots are representative exam- 
ples that illustrate the removal of undesirable high- 
frequency actuator commands; they do not depict 
the PSD for a control law actually tested. The con- 
trol surface rms rate in deg/sec can be calculated as 
the square root of the integral with respect to the 
frequency of the PSD. The design limit for the total 
control surface activity was chosen to be 75 deg/sec, 
as discussed earlier. 



Figure 12. Predicted, control-surface- rate power spectral den- 
sity caused by unit rms gust; closed-loop, symmetric, 
300-psf condition. 


A peak can be seen in the control surface activity 
at a frequency of approximately 11.5 Hz. This 


peak represents the activity required to suppress 
the critical, stabilized, flutter mode as it is excited 
by turbulence. The figure also shows significant 
undesirable control surface activity in the 25-Hz to 
40-Hz frequency range; this activity results from 
turbulence excitation of high-frequency structural 
modes. An analog band-rejection filter was used to 
reduce control surface activity in this region. This 
filter consisted of three fairly broad second-order 
notches with center frequencies at 32 Hz, 40 Hz, 
and 49 Hz. The form of the individual notch elements 
making up the band-rejection filter is shown below: 

BRF = N X N 2 N 3 (8) 

where 

N , = s 2 + %NjUrij S + u%j 
1 - s2 + 

Numerical values for the notch filter parameters 
are shown in the appendix. The band-rejection filter 
was used instead of a low-pass filter to keep the 
resulting lag at the flutter frequency to a minimum, 
while still achieving the desired attenuation. The 
lag at 11.5 Hz, which was due to the filter, was 
approximately 28°. Both table II and figure 12 show 
that the band-rejection filter attenuates a significant 
portion of the undesired high-frequency commands 
to the control surfaces in the 30-Hz to 40-Hz range. 
Prior to the use of the band- rejection filter, the 
rms control rate constraints were violated, and they 
actively contributed to the cost function evaluation 
used in the optimization. After the band-rejection 
filter was designed and incorporated into the design 
plant, the rms control rate constraints no longer 
contributed to the cost function. 

Table II. Predicted Control Rate Reduction With 
Band-Rejection Filter for Modal Rate 
Feedback Design 

[300-psf condition] 



Without filter, 

With filter, 

rms rates 

deg/sec 

deg/ sec 

TEO 

138 

52 

TEI 

11 

4 


See the appendix for the values chosen for pa- 
rameters in the dynamic compensation, the distribu- 
tion matrices ultimately selected, and the blending 
matrices found by the optimization procedure. 
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Predicted ‘performance . Table III shows the 
predicted performance for the modal rate feedback 
controller which resulted from linear analysis and 
which was substantially confirmed by nonlinear batch 
simulation. Both the symmetric and the anti- 
symmetric control laws were predicted to stabilize 
the closed-loop system over the TDT test path to 
the tunnel limit. The gain and phase margins are 
shown, with respect to errors at the single-strand 
point, at the common evaluation point of 300 psf. 
These gain and phase margins were predicted by lin- 
ear analysis to meet the stated requirements through- 
out the wind-tunnel test envelope. The gain mar- 
gins were verified in simulation, at selected dynamic 
pressures, by varying symmetric and antisymmetric 
system gains individually until simulation time histo- 
ries showed divergence. These gain margins obtained 
from simulation were comparable to those obtained 
through linear analysis. The phase margins were not 
verified through batch simulation. 


robustness and performance characteristics. The 
control law was then discretized for implementation. 

Figure 13 shows the analog design plant that con- 
tained the basic state-space model augmented with 
antialiasing filters and a first-order Fade approxima- 
tion for the digital controller 1 computational time 
step delay. The 25-Hz first-order antialiasing filters 
were selected for analysis and design. The TEO and 
LEO control surfaces were used for control input, and 
their collocated accelerometers were used as sensors 
for feedback. Washout filters were added to the con- 
trol law after completion of the modified LQC de- 
sign process. Analog notch filters were also added 
after the design to improve high-frequency robust- 
ness characteristics above 30 Hz. Numerical values 
for controller parameters are presented in the appen- 
dix. In equation form, that part of the control law 
which was implemented digitally is 

6 r = K ; ,(*)z 


Table III. Predicted Performance for Modal Rate 
Feedback Design 


.s 

.S + a 



F) _1 G + e]z 


( 9 ) 


Degrees 
of freedom 

(/max j 

psf 

Margins 
at 300 psf 

rms cont rol 
activity 
at. 300 psf 
(percent of 
max allowed) 

Gain, 

(IB 

Phase, 

<)og 

TEO, 

percent 

TEI, 

percent 

Symmetric 

>325 

±0 

±34 







67 

25 

Antisymmetric 

>325 

±12 

±49 
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The predicted rrris control surface rate was deter- 
mined by using the batch simulation with simulta- 
neous symmetric and antisymmetric turbulence ex- 
citation. The simulation indicated that the specified 
rms control surface rate limit was not exceeded for ei- 
ther pair of control surfaces. Significant activity was 
commanded for both the TEO and the TEI surfaces, 
although the TEO surfaces dominated. 

Modified LQG Design (Design Number 3) 

Overview. Initial symmetric and antisymmetric 
control laws were designed using a modified LQG 
procedure. The philosophy behind this control law 
design was to obtain a minimum energy, full-order, 
optimal controller consisting of a linear quadratic 
regulator and a model-based Kalman state estimator 
for output feedback and then to reduce its order 
without significant loss of the full-order controller 


Design steps. The LQG design and order re- 
duction was arrived at through the iterative proce- 
dure shown in figure 14. Two design points (dynamic 
pressures of 300 psf and 350 psf ) were chosen when' 
the design plant was unstable. The full-state feed- 
back, optimal regulator for each symmetry and de- 
sign point was designed with a zero weighting matrix 
for the states and an identity weighting matrix for 
the controls. For the closed-loop system, this regu- 
lator reflects the unstable plant characteristic roots 
into the left-half plane, while all other roots remain 
unchanged; this represents the minimal control en- 
ergy solution for stabilizing the plant (ref. 56). The 
model-based Kalman state estimator was designed 
with a diagonal fictitious input noise intensity ma- 
trix with elements of 0.000001 rad 2 , a gust input 
noise intensity of 1/144 (ft/sec) 2 , and a diagonal 
measurement noise intensity matrix with elements of 
1/144 (ft/sec 2 ) 2 . 
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Design plant P 3 (s) 



Compensator dynamics 

Figure 13. Block diagram of modified LQG control law. 



Figure 14. Steps in modified LQG design process. 


The full-state feedback regulator was combined 
with the state estimator to generate a full-order com- 
pensator that used only sensor feedback with no 
direct knowledge of the states of the plant. The 


resulting full-order controller required order reduc- 
tion for implementation. The full-order LQG con- 
trol law was reduced through a process of balanced 
realization and modal truncation, based in part upon 
the evaluation of modal residues (ref. 40). Controller 
poles above 25 Hz were removed because they had 
little effect on the control of flutter at 11.5 Hz. A 
lOth-order control law was chosen since its robustness 
and performance characteristics were close to those of 
the full-order LQG control law. The numerical values 
of the parameters in the continuous symmetric and 
antisymmetric control laws, designed for the 300-psf 
and 350-psf points, are shown in the appendix. 
First-order washout filters were added to attenuate 
the response to bias errors, which increased the con- 
troller order to 12. Additional singular value analysis 
revealed the need for stability margin improvement 
in the 32-Hz frequency range. Because the open-loop 
plant poles in this region are stable, signal attenua- 
tion by means of an analog notch filter (see eq. (8)) 
with a center frequency of 32 Hz was added to pre- 
vent the compensator from driving the modes to be 
unstable in this region. 

The final step was to discretize the continuous 
control law at a sample rate of 200 samples per 
second using the Tustin transformation. With this 
transformation, the controller transfer matrix in the 
analog domain and that implemented in the digital 
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Figure IT). Predicted dynamic pressure root loci; symmetric condition; range is from 0 psf to 350 psf at. 50-psf increments. 


domain were virtually identical at and below the 
flutter frequency except for the 1.5 time step de- 
lay present in the digital implementation. The 
Pade approximation to the 1 time step delay that 
was included in the design plant (see fig. 13) re- 
sulted in a 1 time step lead in the continuous con- 
troller that counterbalanced 1 time step of the digital 
implementation delay. 

Poles as function of q. Figure 15 shows 
the predicted plant open-loop poles and predicted 
fixed-gain, closed-loop roots as functions of dynamic 
pressure for the symmetric degrees of freedom. The 
reduced-order 300-psf control law with washouts and 
32-Hz notches was used. To simplify the figure, 
compensator poles and zeros are not shown. Com- 
pensator poles do not change with dynamic pres- 
sure when the feedback loops are open. When feed- 
back loops are closed, the compensator poles interact 
with those of the plant and then change with dy- 
namic pressure. However, for this design, the open- 
arid closed-loop compensator poles are stable for the 
dynamic pressure range shown. 

The solid lines in figure 15 indicate the paths of 
the open-loop poles, and the dashed lines indicate the 
paths of the closed-loop roots. The crossing point 
where the 11.5-Hz flutter mode becomes unstable is 
identified in the figure as 248 psf for the design model 
of the symmetric plant with no compensation and as 


350 psf for the symmetric plant with compensation. 
The corresponding values for the antisymmetric flut- 
ter mode were 252 psf and 325 psf, respectively. (The 
root loci are not shown.) 

Predicted performance . The modified LQG 
control law based on the 300-psf design point was 
predicted through analysis and simulation to provide 
closed-loop stability to the limit of the wind-tunnel 
operating range. The gain and phase margins shown 
in table IV represent guaranteed minimum margins 
for simultaneous variations on multiple channels 
(ref. 57). These margins can be conservative if they 
represent an unlikely combination of variations. The 
margins shown here do not meet the requirement s for 
SISO gain and phase margins. However, because of 
their potentially conservative nature, these margins 
were judged to be sufficient for testing the control 
law. 

The closed-loop rms control surface rates in the 
presence of random gust excitation are within the 
specified limits. The percent of maximum allowed 
control surface activity for each pair of surfaces in- 
dicates that both the TEO and LEO control surface 
pairs are used to a significant extent, but the TEO 
surfaces are dominant. The rms control surface rates 
were generated with both symmetries active by us- 
ing the batch simulation. Separate linear analyses 
of rms control rate activity were performed for each 
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Table IV. Predicted Performance for Modified LQG Design 
[Laws based on 300-psf design point] 


Degrees 
of freedom 

^max i 

psf 

Margins 
at 300 psf 
(at plant 
input)" 

rms control 
activity 
at 300 psf 
(percent of 
max allowed) 

Gain, 

dB 

Phase, 

deg 

TEO, 

percent 

LEO. 

percent 

Symmetric 

350 

±3 

±18 







62 

26 

Antisymmetric 

325 

±4 

±20 




"Equivalent multivariable margins for simultaneous inde- 
pendent changes on all channels. 


symmetry for the 300-psf controller. These analy- 
ses showed that, for a unit gust intensity, the anti- 
symmetric control activity was only approximately 
25 percent as large as the symmetric control activity. 

Wind-Tunnel Test Results 

Measured Versus Predicted AFW 

Wind-Tunnel Model Characteristics 

A variety of information was collected during the 
test. Runs that were made early in the entry es- 
tablished that the uncontrolled flutter dynamic pres- 
sure of the decoupled tip ballast store configuration 
was sufficiently high to provide the desired flutter- 
stopper capability. This capability was achieved, 
however, only after adjusting the stiffness of the tor- 
sional spring that is activated in the decoupled con- 
figuration to avoid adverse coupling with model elas- 
tic modes of approximately 6 Hz. The stiffness was 
adjusted to reduce the frequency of a tip ballast store 
mode from 6 Hz to 4.5 Hz. The possibility of this 
requirement had been anticipated, and the stiffness 
adjustment capability had been built into the design. 

Flutter clearance runs were made with the tip 
ballast stores in the coupled configuration to estab- 
lish the uncontrolled flutter boundary for the con- 
figuration that was to be tested in a closed-loop 
manner. Differences were observed between pre- 
dicted and actual flutter dynamic pressures. The dy- 
namic pressure for antisymmetric flutter was found 
to be lower than that predicted by approximately 
30 psf or 13 percent based upon a large, primarily 
antisymmetric response encountered at a dynamic 
pressure of about 220 psf. Analysis had indicated 
that symmetric flutter would occur first, at approxi- 
mately 248 psf, with antisymmetric flutter occurring 
at about 252 psf. 


Subcritical open-loop runs were made to obtain 
estimates of plant frequency response functions. In 
these runs, excitation was generated within the dig- 
ital computer and sent out to an actuator pair as 
either a symmetric or an antisymmetric command. 
A sinusoidal excitation was input with the frequency 
varying logarithmically with time over a 150-sec pe- 
riod from 4 Hz to 35 Hz. The input and outputs 
of interest were recorded digitally and processed, in 
near real time, using fast Fourier transform tech- 
niques to obtain plant estimates. A key difference 
between predicted and actual characteristics was ev- 
ident. For both symmetries, the frequencies at which 
dominant frequency response peaks occurred were 
somewhat lower than predicted. Figure 16(a), which 
shows the symmetric case at a dynamic pressure of 
175 psf, contains a comparison of predicted and mea- 
sured Bode plots of ^tip/^TECV This sensor/control 
combination exhibits similarities and differences that 
are typical of what has been seen with other combi- 
nations. Both curves contain the effect of the anti- 
aliasing filters and a 0.5 time step delay. (The 1 time 
step delay associated with controller output com- 
mands is not present here.) The predicted and mea- 
sured curves exhibit the same trend over the fre- 
quency range shown, and they show good agreement 
in peak magnitudes. However, a shift of approxi- 
mately 1 Hz is evident in the frequency at which 
the peak occurs for the lower frequency of the two 
interacting modes. The analytical model had been 
adjusted so that the frequencies at zero dynamic 
pressure matched the frequencies measured during a 
ground vibration test. The differences between pre- 
dicted and measured frequencies must, therefore, be 
related to aerodynamic effects that could arise from 
a number of sources, including errors in the predicted 
mode shapes. 

The phase characteristics of the response shown 
in figure 16(a) indicate a frequency shift consistent 
with the frequency shift for the peak magnitude. 
(The predicted phase at the predicted frequency of 
peak response is in close agreement with the actual 
phase at the actual frequency of peak response.) 
Figure 16(b), which presents a polar plot of the 
curves of figure 16(a), illustrates the correspondence 
between the magnitude and phase frequency shifts. 
The phase correspondence between the measured and 
predicted curves is good, particularly in the vicinity 
of the peak response. 

If a control law has dynamics in the flutter 
frequency range, the shifted plant dynamics can 
introduce potentially large phase shifts in the fre- 
quency response of the loop transfer function. There- 
fore, control law designers should be aware of the 
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Measured Predicted 




(b) Polar plot. 

Figure 1 6. Predicted and measured open-loop frequency re- 
sponse for 2 ] ip resulting from <^teC) ( • symmetric. 175-psf 
condition. (Predicted response includes effects of 25-Hz 
antialiasing filter and 0.5 time step delay.) 


sensitivity of their designs to changes in the frequen- 
cies of critical modes and should not rely strictly on 
the adequacy of typical phase margin requirements. 


Because none of the control laws was scheduled 
with dynamic pressure, it is more significant to com- 
pare the difference between measured and predicted 
frequencies at dynamic pressures that are the same 
percent away from the corresponding flutter dynamic 
pressures than it is to compare the difference be- 
tween measured and predicted frequencies at a given 
dynamic pressure. The relevant frequency shift at 
flutter was approximately 2 Hz. 


Traditional Pole/Zero Loci Controller 

The traditional polc/zcro design demonstrated 
closed-loop stability up to a dynamic pressure of 
about 272 psf. This test represented an increase 
of approximately 24 percent relative to the observed 
open-loop flutter boundary. Flutter was suppressed 
simultaneously in both symmetries. The controller 
stabilized the model at the 272-psf condition as in- 
dicated by the decay of bursts of turbulence-induced 
system response. The wind tunnel was operated at 
this condition for several minutes while time histories 
for loads and for commanded control deflections re- 
sulting from tunnel turbulence were being recorded 
for rms analysis. The rms control rates were only 
about 25 deg/scc, which is one-third of the accept- 
able maximum. The wind-tunnel safety system was 
activated automatically after the model responded 
to a burst of larger amplitude turbulence and the 
structural loads exceeded preset limits. Because the 
control law stabilized the system and was also able 
to limit the amplitude of the flutter mode for lower 
turbulence levels with significant, reserve rate capabil- 
ity, it is speculated that increasing the feedback gain 
would have kept the structural loads caused by tur- 
bulence to be within the prescribed limits, at least in 
the flutter frequency range; however, a reduced high- 
frequency gain margin and, possibly, smaller phase 
margins would have accompanied the feedback gain 
increase. 

Figure 17 presents singular value assessments of 
the antisymmetric control law obtained using the 
controller performance evaluation (CPE) software 
(ref. 58) and experimental data. This control law 
provides two inputs to the plant and utilizes one sen- 
sor output. Figures 17(a) and 17(b) present singular 
values of the return difference matrices at the plant 
input and plant output, respectively. Large values (a 
value of 1 is large) for both sets of minimum singu- 
lar values would indicate that the closed-loop system 
stability characteristics are tolerant to unstructured 
uncertainties at the plant input and output. 

The relatively small minimum singular values at 
the plant input which are seen in figure 17(a) near 
7 Hz would result from errors in the worst possi- 
ble direction and arc not necessarily likely or even 
physically realizable. Consequently, use of this fig- 
ure to assess tolerance to errors may be quite con- 
servative. Error sources that could occur include er- 
rors in individual control surface aerodynamic effec- 
tivenesses in the 7-Hz frequency region, in control 
surface mass coupling terms, and in gain and phase 
of commands from the controller to the individual 
actuators. The successful closed-loop tests demon- 
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(a) Return difference at plant input; max is 9.69; min is 0.29 
at 7.1 Hz. 



(h) Return difference at plant output; max is 4.71; min is 0.51 
at 6.7 Hz. 


Figure 17. Experimentally determined singular values for tra- 
ditional pole/zero loci control; antisymmetric, 240-psf 
condition. 

strated that the small but conservative assessment of 
tolercncc to uncertainty of figure 17(a) was sufficient. 

Only one singular value curve exists in fig- 
ure 17(b) for the scalar plant output, and its mag- 
nitude corresponds to the distance of the single- 
strand Nyquist plot from the singular point. The 
singular values of figure 17(b) also correspond to 
the distance of the multi-input /single-output (MISO) 
Nyquist plot from the singular point. Figure 17(b) 
does not address distinct error sources in the two in- 
put channels, and, therefore, it is a potentially non- 
conservative assessment of robustness. Nevertheless, 
the test results demonstrated that the use of SISO ro- 
bustness criteria to assess MISO robustness was ade- 
quate, in this case, to obtain a successful design. The 


global minimum near 7 Hz in figure 17 was associated 
with sensitivity to uncertainty in phase lead. 

Modal Rate Feedback Controller 

The modal rate feedback controller that was de- 
signed based upon the predicted AFW wind-tunnel 
model characteristics was shown, by experimentally 
derived open-loop CPE, to be destabilizing at the two 
highest subcritical dynamic pressures tested (125 psf 
and 175 psf). The primary cause for these insta- 
bilities was believed to be undue sensitivity of the 
controller dynamics to frequency shifts of the critical 
structural modes. 

Because the design method can readily employ 
experimentally derived frequency responses as inputs 
to the optimization, the frequencies of the controller 
dynamics were shifted to match the observed shift, 
and the blending matrix was reoptimized using trans- 
fer matrix estimates based upon data collected at 
125 psf and 175 psf. No constraints were placed 
upon rms controller rate requirements because no ex- 
perimentally derived frequency responses existed for 
outputs due to gust inputs. The numerical values of 
the parameters in the reoptimized control laws are 
shown in the appendix. Subsequent open-loop CPE 
and closed-loop testing with the redesigned controller 
showed that the system performance at 125 psf and 
175 psf was in agreement with what had been pre- 
dicted using the earlier experimental data. However, 
the controller destabilized the system at 185 psf with 
the instability occurring at a frequency of approx- 
imately 7 Hz. The source of the deficiency is dis- 
cussed subsequently in the section entitled “Posttest 
Analysis.” 

Modified LQG Controller 

The modified LQG controller designed for 300 psf 
did not significantly change the closed-loop flutter 
dynamic pressure relative to the observed open-loop 
dynamic pressure. With this controller operating, 
antisymmetric flutter, at approximately 9.5 Hz, was 
encountered near the observed open-loop boundary 
of 220 psf. Data from the CPE analysis indicated a 
much lower antisymmetric component of control sur- 
face activity than predicted, which raised the possi- 
bility that the gain for the antisymmetric control law 
was too low. This aspect is discussed further in the 
section entitled “Posttest Analysis.” 

The controller designed for 350 psf was also 
tested. This 350- psf controller was very similar to 
the 300-psf controller, but it operated at higher gain 
levels. With the 350-psf controller operating, anti- 
symmetric control surface activity was significantly 
higher, but the symmetric control law gain was too 
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high in the 21-Hz region. The closed-loop system 
was driven symmetrically unstable at 175 psf at a 
frequency of approximately 21 Hz. Additional dis- 
cussion of the sources of these difficulties is presented 
in the section entitled “Posttest Analysis.' 1 

Combined Performance 

Figure 18 shows the maximum closed-loop sta- 
ble dynamic pressure achieved experimentally by 
each control law. Figure 19 depicts closed- loop con- 
trol surface rms rates as functions of dynamic pres- 
sure. The rms rates were estimated bv differenti- 
ating commanded deflections because rates were not 


O Traditional pole/zero loci 



Figure 18. Maximum dynamic pressure obtained during 
closed-loop testing. 



Figure 19. Experimentally determined root- mean-square (rms) 
control rates (design limit was 75 deg/sec rms). 


commanded directly. The three plots represent the 
measured response for the three control laws during 
testing. Because the TEO control surfaces were 
dominant for each control law, rates are shown for 
the TEO surfaces. 

All three control laws command similar levels 
of rms control surface activity at common dynamic 
pressures tested. This similarity reflects the fact that 
all three were designed with the same turbulence 
model and the same design limits. None of the con- 
trol laws had difficulty staying within the design limit 
of 75° /sec rms. In fact, the peak measured rms rate is 
only approximately one-third of the limit for the tra- 
ditional pole /zero design at 272 psf. The difference 
between predicted (see tables I, III, and IV) and mea- 
sured rms control rate requirements may be because 
of errors in assumed turbulence levels and models, in 
control effectivenesses, and in modal damping. The 
low control surface rms rates indicate that additional 
control power is available for improved performance. 

Posttest Analysis 

Modal Rate Feedback Design 

The test results clearly indicated that the pretest 
design and analysis associated with robustness as- 
sessments with respect to errors at the single-strand 
point (SISO) were insufficient. The basic problem 
was that, although there was an SISO point in the 
feedback loop, multiple sources existed for gain and 
phase errors, that is, four pairs of sensors and two 
pairs of controls. The optimization procedure did 
not consider control law performance sensitivities to 
multiple, independent errors. Singular value anal- 
ysis of errors at the plant input and output points 
(fig. 9 shows the loop breaking points as dotted ovals) 
provides a means, although potentially quite conser- 
vative, to assess this sensitivity. Figure 20 graph- 
ically illustrates the deficiency of the design. This 
figure was generated using experimentally derived 
plant transfer matrix elements and the control law 
developed using the experimental data. This figure 
presents logarithmic plots of the minimum singular 
values of the return difference matrices at the single- 
strand point, the plant input point, and the plant 
output point for both symmetries. Large minimum 
singular values are desirable for robustness to error. 
Extremely small minimum singular values are evi- 
dent for the return difference matrix at the plant 
output in a frequency range near 7 Hz. 

The small minimum singular values were a result 
of the choice of objective function for controller pa- 
rameter optimization. Consider the objective func- 
tion (see eq. (7)) in the typical condition in which the 
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(a) Return differences at symmetric condition. 



(b) Return differences at antisymmetric condition. 

Figure 20. Experimentally determined minimum singular val- 
ues for modal rate feedback control; 175-psf condition; law 
is derived from experimental data. 

control rms constraints were inactive. In this condi- 
tion, the optimization procedure was attempting to 
increase the nonconservativc, upper bound SISO sin- 
gular value abs{l +Tgs(.s)} (ref. 35), and it was al- 
lowed to do so at the expense of the MIMO minimum 
singular values at the input a\I 4- K‘ 2 (s)P 2 (s)] and 
the output <t[1 + P 2 (s)K 2 (s)]. Extreme sensitivity of 
the closed-loop system to errors at the plant output 
resulted and led to destabilization by the controller. 
This result emphasizes the importance of properly 
capturing all the critical design trade-offs in either 
the objective function or in the design constraints. 

Modified LQG Design 

The reduced-order controller developed using this 
design process and based upon a design point at 
300 psf did not raise the closed-loop flutter dynamic 
pressure. The basic deficiency in the controller is 


shown in figures 21 and 22. Figure 21 shows two plots 
of the determinant of the matrix [I + P 3 (s)K 3 (s)], 
which is the MIMO Nyquist plot with zero as the 
critical point. The dashed curve shows predicted val- 
ues, and the solid curve shows values derived from 
the experimental CPE analysis; both curves repre- 
sent the antisymmetric condition at a dynamic pres- 
sure of 200 psf. The control law employed was the 
one based on the 300-psf design point. The small- 
ness of the experimental curve relative to the ana- 
lytical one is surprising. The test point of 200 psf 
is only 20 psf from the observed open-loop flutter 
dynamic pressure, whereas, for the analytical predic- 
tions, 200 psf is 52 psf from open-loop flutter. One 
would, therefore, expect the plant response to in- 
crease in amplitude much more in the test than in 
the analytical predictions. The indication is that the 
controller gain is too low. Figure 22 presents fre- 
quency responses for each channel of the continuous 
form of the reduced-order controller for each sym- 
metry. (The washout filters are not included.) The 
antisymmetric control law has low gain in all chan- 
nels in the 9-Hz to 10-Hz flutter frequency range as 
compared with the symmetric law. This difference in 
gain level is in contrast to the pole/zero loci law for 
which the dominant channel (^TEO. /^TIp) was 
same for both symmetries. The gain in the antisym- 
metric channels was not sufficient to suppress flut- 
ter; however, the size of the antisymmetric controller 
magnitude relative to the symmetric controller mag- 
nitude does not fully explain the drastic difference 
between the predicted and measured MIMO Nyquist 
curves in figure 21. 



Real 

Figure 21. Predicted and measured multi-input/multi-output 
Nyquist plot for LQG law; antisymmetric, 200-psf condi- 
tion; law is based on 300-psf design point. 
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Frequency, Hz 

(b) Antisymmetric. 

Figure' 22. Frequency response for modified LQG controller 
designed at 300 psf. 

The frequency responses for the controller design 
based upon the 350-psf point are shown in figure 23 
and exhibit similar but higher gain characteristics as 
compared with the 300-psf design. Note the large 
amplitudes in the region of 21 Hz, particularly for 
^TEO< / ^TEO f° r the symmetric law. This peak is 
present as a result of the recovery process and the ex- 
istence of a lightly damped nonminimum phase trans- 
mission zero in the symmetric state-space model of 
the plant. The lowly damped controller pole is near 
the mirror image (with respect to the imaginary axis) 
of the nonminimum phase plant transmission zero. 


Such a factor makes the design susceptible to uncer- 
tainty in the plant characteristics near 21 Hz, which 
is a frequency approximately double the frequency to 
be controlled. This sensitivity is shown in figure 24, 
which presents predicted minimum singular values 
associated with the additive plant error at 200 psf 
for the symmetric degrees of freedom; the control law 
designed for the higher gain. 350-psf design point, 
was used. The figure reveals a controller stability 
robustness sensitivity to error not considered in the 
design process. An undesirably low minimum value 
occurs at a frequency near 21 Hz. 

S TEO c S LEO c ^TEO c 5 LEO c 

z LEO z LEO z TEO z TEO 



(a) Symmetric. 


5 leo c s teo c s leo c ^teo c 

z TEO z TEO z LEO z LEO 



(b) Antisymmetric. 

Figure 23. Frequency responses for modified LQG controller 
designed at 350 psf. 
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Figure 24 . l > redicted inininmm singular value for modified 
LQG control; symmetric, 200-psf condition; law is based 
on 350-psf design point. 


The high controller gain near 21 Hz, coupled 
with error between the predicted and actual plant, 
destabilized the closed-loop system during the test, 
near 175 psf. It is conjectured that better per- 
formance would have been achieved with the anti- 
symmetric law designed based upon the 350-psf point 
and the symmetric law designed based upon the 
300- psf point. 

Concluding Remarks 

Three flutter suppression control laws were de- 
signed in the continuous domain for the Rockwell In- 
ternational Corporation active flexible wing (AFW) 
wind-tunnel model. The control laws were imple- 
mented digitally and tested subsonically in the Lang- 
ley Transonic Dynamics Tunnel at Langley Research 
Center. All three control laws were predicted to meet 
the objective of significantly raising the flutter onset 
dynamic pressure while not violating control surface 
rate and displacement limits. 

Open-loop wind-tunnel testing exposed differ- 
ences between predicted and actual AFW wind- 
tunnel model characteristics, particularly in the fre- 
quency of the flutter mode. However, the analytical 
model correctly captured the essential character of 
the flutter mechanism. 

Only the traditional polc/zero loci design was 
sufficiently robust to the model errors to raise the 
closed-loop flutter dynamic pressure. With the 
pole/zero loci design, simultaneous suppression of 
symmetric and antisymmetric flutter was success- 
fully demonstrated to a dynamic pressure 24 percent 
above the open-loop boundary. At this condition, the 


controller still provided stability as demonstrated by 
its succuessful damping of response resulting from 
typical bursts of turbulence over a period of 2 min- 
utes; subsequently, however, a burst of larger am- 
plitude turbulence caused torsional loads to exceed 
preset safety limits, at which point testing of this 
law was terminated. 

The multi-input/multi-output modal rate feed- 
back controller design process did not incorporate 
into the optimization procedure constraints upon 
critical sensitivities to errors at the plant output and 
to frequency shifts in the plant dynamics. This de- 
sign led to inadequate robustness to the modeling 
errors that were encountered and closed-loop insta- 
bility at a lower dynamic pressure than that of the 
observed open- loop condition. 

The modified linear quadratic Gaussian controller 
design with subsequent controller order reduction 
was also sensitive to design model errors and did 
not have adequate guaranteed gain and phase mar- 
gins. The sensitivity and robustness characteristics 
were strongly influenced by the presence of a lightly 
damped nonminimum phase transmission zero in the 
design model of the plant. As implemented, the de- 
sign process placed a controller pole near the mirror 
image location, with respect to the imaginary axis, 
of the nonminimum phase zero. This placement not 
only resulted in undesirably high gain and sensitivity 
to plant error at a frequency that was double that 
of the mode to be controlled but also constrained 
the gain in the flutter frequency region. Closed- 
loop flutter was encountered during the test, near 
the open- loop flutter dynamic pressure at the uncon- 
trolled flutter frequency for a low-gain controller arid 
below the open-loop flutter dynamic pressure at the 
frequency of the design model nonminimum phase 
transmission zero for a high-gain controller. 

The November 1989 test provided data for as- 
sessing the fidelity of the analytical models of the 
AFW wind-tunnel model and for evaluating the ro- 
bustness of the control laws to real-world implemen- 
tation considerations. The lessons that were learned 
were applied in a subsequent effort in which four 
separate flutter suppression control laws were suc- 
cessfully tested not only in steady flight but also 
while performing aggressive, actively controlled, roll 
maneuvers. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
July 21, 1992 


28 



Appendix 

Numerical Definition of Controller Parameters 

The numerical values are specified here for the parameters in the continuous domain representation of the 
control laws that have been developed using each of the three design approaches. These definitions correspond 
to control laws tested in the Langley Transonic Dynamics Tunnel in November 1989. In each case, a digital 
implementation was carried out, prior to testing, using a Tustin transformation with no prewarping and a 
sample rate of 200 Hz. 


Traditional Pole/Zero Loci Control Law Definition 

The parameters (see fig. 6 for their significance) have the same values for both symmetries. For the 
symmetric case, the switch is open in the channel commanding the trailing-edge inboard actuator pair. 


a = 5.000 rad/sec (0.7958 Hz) 


Cjv = 0.4706 

UJ N = 85.00 rad/sec (13.53 Hz) 

Co = 0.09950 

uj d = 70.35 rad/sec (11.20 Hz) 

k\ = 0.4871 deg/g (streamwise) 

d = -0.2500 


This design did not consider the effect of computational delay. Consequently, after applying the Tustin 
transformation, the “buy-back” procedure (see eqs. (2) and (3)) was employed to approximately counter the 
effect of the delay. 


Modal Rate Feedback Control Law Definition 


Figure 9 shows the significance of the parameters. 


AAF 

BRF 

N i 


N 2 

n 3 

T 

INTI 

INT2 

WOF 

DNOT 


157 

,s + 157 
N 1 N 2 N 3 

•S 2 + 2(0.08) (200)* + 200 2 

s 2 + 2(0.32)(200).s + 200 2 

,s 2 + 2(0.16)(250).s + 250 2 

,s 2 + 2(0.48)(250).s + 250 2 

,s 2 + 2(0.12)(310 )j» + ,'ilO 2 

T 2 + 2(0.32)(310)» + 310 2 

1.5T = 0.0075 sec 
1 

s -f «| 

1 

,s + «2 
,s 

.s -f a 

4 + ijjf, 

.s 2 4 2( t i)n) n s 4 


(25.0 Hz) 


(31.8 Hz) 
(39.8 Hz) 
(49.3 Hz) 
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i 

s 

N 

l 

*r 

= B ln 

^TEI 




2TEO 



- ^TIP - 


^TEI { . 

. S TFjO c . 


= VisU 


Analytically derived control law . 

Symmetric: 


a i 

= 

4 

rad/sec 


= 

30 

rad /sec 


= 

40 

rad/sec 


= 

36 

rad /sec 


(0.6 Hz) 
(4.8 Hz) 
(6.4 Hz) 
(5.7 Hz) 


-0.0482 0.0780] * 70 

-0.5709 0.4476] * 70 2 


Cv = 

C D = 


B 


LN = 


D 


IS = 


0.03 

0.30 

"[-0.3857 0.1187 
[-0.5276 1.0000 

"0.4500" 

1.0000 


k ,2 - -9.9251 deg/ (g-scc) (degrees are streamwise) 


Antisymmetric: 


a ] = 

a 2 = 
a = 
Vn = 


Cv = 

Cd = 

Bln = 

Dis = 


k 2 = 


6 

rad/ sec 

(1.0 

Hz) 


56 

rad / sec 

(8.9 

Hz) 


62 

rad /sec 

(9.9 

Hz) 


115 

rad/sec 

(18.; 

i Hz) 


0.03 





0.10 





'[-0.5480 0.3341 - 

•0.1529 

0.2053] 

[-0.4165 1.0000 - 

0.7755 

0.5353] 


0.9000" 

1.0000 

—4.1000 dcg/(g-sec) (degrees are streamwise) 
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Experimentally derived control law. 
Symmetric: 


a \ - 

4 rad/sec 

(0.6 Hz) 


«2 = 

1 7 rad / sec 

(2.7 Hz) 


a — 

25 rad/sec 

(4.0 Hz) 


— 

33 rad /sec 

(5.3 Hz) 


Cat = 

0.06 




II 

0.15 




b ln = 

[1.0000 

[0.3196 

0.0232 

0.2871 

-0.4618 

-0.0060 

-0.1479] 

-0.2931] 

D,s = 

"0.4500" 

1.0000 





* 70 

* 70 2 


A ,’2 = 2.5000 dog/ (g-sec) (degrees are streamwise) 


Antisymmetric: 




a { = 

6 rad /sec 

(1.0 Hz) 


«2 = 

38 rad / sec 

(6.0 Hz) 


a — 

42 rad /sec 

(6.7 Hz) 


— 

103 rad /sec 

(16.4 Hz) 


Cv = 

0.03 




C D = 

0.30 




B ln = 

'[ 0.1578 -0.0867 -0.1723 

[-0.3481 1.0000 -0.1156 

0.0677] * 70 

-0.1151] * 70 2 

Bis = 

"0.9000" 

1.0000 




k ‘2 — 4.3000 deg/(g-sec^ 

) (degrees 

are streamwise) 


Modified LQG Control Laws 

These matrices define the continuous controllers designed by the modified linear quadratic Gaussian (LQG) 
procedure followed by controller order reduction. Controller results are shown for design points at 300 psf and 
350 psf. The design point upon which the controller is based is included as a subscript on the state matrices. 
Likewise, the subscript S refers to symmetric and A refers to antisymmetric. These matrices were designed 
with a Pade approximation of a 0.005-sec time delay included as part of the design plant (fig. 13). Figure 13 
also shows elements added after the LQG design. These elements are analog notch filters ( N\ as defined 
in the section entitled “Modal Rate Feedback Control Law Definition' 1 in this appendix) and washout filters 
(a = 6 rad/sec). The washout filters increase the digitally implemented order of the controller portion by 2 
to 12. The inputs to the control laws are accelerations in gravitational units, and the outputs are commanded 
control deflections in degrees streamwise. 
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F3005 = 


■ -10.2606 

11.4264 

0 

0 

0 

0 

0 

0 

0 

0 

-11.4264 

-10.2606 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.8167 

36.8465 

0 

0 

0 

0 

0 

0 

0 

0 

-36.8465 

-1.8167 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-29.4794 

51.0219 

0 

0 

0 

0 

0 

0 

0 

0 

-51.0219 

-29.4794 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-25.7801 

113.8901 

0 

0 

0 

0 

0 

0 

0 

0 

-113.8901 

-25.7801 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-6.0074 

132.0979 

0 

0 

0 

0 

0 

0 

0 

0 

-132.0979 

-6.0074. 

G3005 = 










■ 0.1044 

0.1176 “ 









-3.5706 

-5.6319 









-1.4527 

0.7643 









0.6457 

-0.2196 









5.6922 

-4.3832 









-4.6252 

2.4555 









3.9561 

-3.8048 









-0.1957 

0.4442 









0.7371 

0.5725 









.-1.1023 

-3.1639. 










H3(% = 

-6.9516 -0.6589 

8.7253 -1.5487 


-3.0557 1.3527 

3.0937 -5.5602 


-5.9818 1.0610 

10.7247 -7.8683 


8.6213 -0.2974 

-22.1874 9.6035 


4.2282 1.8143 

-14.2105 -0.9051 


E30O5 = 

'-0.0710 0.0003 
0.2204 0.0479 
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F.'ioo.i = 


--8.3952 

8.9838 

0 

0 

0 

0 

000 

0 

-8.9838 

-8.3952 

0 

0 

0 

0 

000 

0 

0 

0 

-28.5058 

61.1189 

0 

0 

000 

0 

0 

0 

-61.1189 

-28.5058 

0 

0 

000 

0 

0 

0 

0 

0 

-2.8250 

111.6667 

0 0 0 

0 

0 

0 

0 

0 

-111.6667 

-2.8250 

000 

0 

0 

0 

0 

0 

0 

0 

-28.0111 122.2674 0 

0 

0 

0 

0 

0 

0 

0 

-122.2674 -28.0111 0 

0 

0 

0 

0 

0 

0 

0 

0 0 -8.2651 164.0179 

0 

0 

0 

0 

0 

0 

0 0 -164.0179 - 

-8.2651 

G3004 = 








- -0.0364 

-0.0583“ 







1.9539 

3.1044 







-2.1703 

1.8668 







1.9796 

-1.0803 







-0.1316 

0.2361 







0.8289 

0.0775 







-2.4537 

2.7932 







-1.5799 

1.2050 







-0.0920 

-0.8088 







_ 4.8785 

-0.6457. 







H 3 oo 4 = 








2.5501 

0.0272 

5.3665 

-1.9025 

-0.9993 

0.1522 - 

-4.2701 0.5809 -1.2457 0.1270' 


-5.9806 

0.5921 

-15.6204 

9.3613 

3.4731 - 

1.1375 

14.5587 -4.3242 4.4800 -1.0030 


= 








[-0.0018 

-0.0152" 

1 







0.0083 0.0617 
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F;i50 s . = 


' -9.5757 

10.8686 

0 

0 

0 

0 

0 

0 0 

0 

-10.8686 

-9.5757 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

-2.2954 

36.7980 

0 

0 

0 

0 0 

0 

0 

0 

-36.7980 

-2.2954 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

-25.8883 

41.6699 

0 

0 0 

0 

0 

0 

0 

0 

-41.6699 

-25.8883 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

-20.9726 

111.5183 0 

0 

0 

0 

0 

0 

0 

0 

-111.5183 

-20.9726 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 -5.2954 132.0236 

0 

0 

0 

0 

0 

0 

0 

0 -132.0236 - 

5.2954 

G;ir,o s = 

' 0.5227 

—0.2915 " 








-4.1499 

-7.3690 








-4.7915 

2.3804 








0.0094 

0.5254 








7.7598 

-5.1526 








-3.8145 

1.7785 








-3.3217 

3.1583 








2.3376 

-2.5772 








2.0936 

1.3330 








. -0.9683 

-8.0868. 








H 35() 6 - = 
'-10.0640 

-1.4205 

-2.6514 

-0.1264 

-5.8140 - 

0.3791 - 

- 9.7835 

1.4382 1.9375 1.6874] 


8.4008 

-2.2960 

2.6435 

-3.5681 

7.2858 - 

-9.7980 

16.8243 -19.1192 -8.5543 -0.9016 J 


E 350 5 = 
f -0.0849 

0.0238] 









0.2998 0.0567 
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F.'irjii, = 


—7.6193 

8.5074 

0 

0 

0 

0 

0 0 0 

0 

-8.5074 

-7.6193 

0 

0 

0 

0 

0 0 0 

0 

0 

0 

-28.6342 

53.7216 

0 

0 

0 0 0 

0 

0 

0 

-53.7216 

-28.6342 

0 

0 

0 0 0 

0 

0 

0 

0 

0 

-2.8832 

111.9853 

0 0 0 

0 

0 

0 

0 

0 

-111.9853 -2.8832 

0 0 0 

0 

0 

0 

0 

0 

0 

0 

-22.3533 120.0340 0 
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0 
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0 

0 

0 

0 

1 

to 
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to 

to 

yC 
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0 

0 
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0 
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0 0 -7.1933 

163.7060 
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0 0 -163.7060 

-7.1933 

G;*r,0 4 = 








'—0.2264 

—0.0472 “ 







3.0473 

5.1245 







-5.8241 

4.5198 







3.7363 

-1.9123 







-0.1919 

0.4602 







1.2983 

-0.2419 







2.9017 

-3.5133 







1.9497 

-1.7576 







-0.1752 

-0.9363 







. 6.8011 

-0.8170. 







H;350^ = 
4.8726 

0.5716 

3.7646 

-0.2633 

-1.2795 

0.1815 5 

.0504 0.9597 -1.4415 -0.12211 


-8.6031 

0.8148 

-8.8879 

7.9763 

4.8913 - 

3.1795 -19.8503 6.3186 6.1731 -1.6116' 


E;350 4 = 








[-0.0106 

-0.0132' 

1 







0.0792 0.1139 
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